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Abstract. Neutralino annihilation in the Galactic halo is the most definite observational 
signature proposed for indirect registration of the SUSY Dark Matter (DM) candidate particles. 
The corresponding annihilation signal (in the form of gamma-rays, positrons and antiprotons) 
may be boosted for one or three orders of magnitude due to the clustering of cold DM particles 
into the small-scale and very dense self-gravitating clumps. We discuss the formation of these 
clumps from the initial density perturbations and their successive fate in the Galactic halo. Only 
a small fraction of these clumps, ~ 0.1 %, in each logarithmic mass interval A log M ~ 1 survives 
the stage of hierarchical clustering. We calculate the probability of surviving the remnants of 
dark matter clumps in the Galaxy by modelling the tidal destruction of the small-scale clumps 
by the Galactic disk and stars. It is demonstrated that a substantial fraction of clump remnants 
may survive through the tidal destruction during the lifetime of the Galaxy. The resulting mass 
spectrum of survived clumps is extended down to the mass of the core of the cosmologically 
produced clumps with a minimal mass. The survived dense remnants of tidally destructed 
clumps provide an amplification (boosting) of the annihilation signal with respect to the diffuse 
DM in the Galactic halo. We describe the anisotropy of clump distribution caused by the tidal 
destruction of clumps in the Galactic disk. 


1. Introduction 

The stable SUSY particle like neutralino is a promising candidate for the enigmatic Dark Matter 
(DM) particle. A primordial power-law spectrum of density fluctuations in the DM ranges from 
the largest scales above the scales of superclusters of galaxies to the smallest sub-stellar scales 
according to prediction of inflation models. This permits to predict the properties of smallest 
DM structures from the known Cosmic Microwave Background (CMB) fluctuations at large 
scales. The most bright indirect signature of the SUSY DM particles is their annihilation, and 
this process could be boosted inside the dense DM clumps in the galactic halo. 

The cosmological formation and evolution of small-scale DM clumps have been studied in 
numerous works [DGlGlElGliElIZlElEl Dm IHJ 1121 H3 III ESI I16| . The minimum mass of 
clumps (the cutoff of the mass spectrum), M m \„ is determined by the collision and collisionless 
damping processes (see, e. g., [5J and references therein). Additionally the cutoff of mass 
spectrum is influenced by the acoustic absorption m at the time of kinetic decoupling of 
the SUSY DM particles [18] and also by the horizon-scale perturbation modes m- The low- 
mass cut-off of the clump mass-spectrum accompanies the process of decoupling. It starts when 
DM particles coupled strongly with surrounding plasma in the growing density fluctuations. 
The smearing of the small-scale fluctuations is due to the collision damping occurring just 


before decoupling, in analogy with the Silk damping [20]. It occurs due to diffusion of DM 
particles from a growing fluctuation, and only the small-scale fluctuations can be destroyed 
by this process. The corresponding diffusive cut-off M is very small. As coupling becomes 
weaker, the larger fluctuations are destroyed and increases. One may expect that the 

largest value of M m [ n is related to a free-streaming regime. However, as recent calculations 
show [|9j, the largest M m - m is related to some friction between DM particles and cosmic plasma 
similar to the Silk damping. The predicted minimal clump masses range from very low values, 
M mir, ~ 10 _12 Mq [21], produced by diffusive escape of DM particles, up to M m - m ~ 10 -4 M Q , 
caused by acoustics oscillations m and quasi-free-streaming with limited friction [S], In the 
case of the Harrison-Zeldovich spectrum of primordial fluctuations with CMB normalization the 
first Earth-mass small-scale DM clumps are formed at redshift z ~ 60 (for 2 a fluctuations) with 
a mean density 7 x 1CD 22 g cm -3 , virial radius 6 x 10 -3 pc and internal velocity dispersion 
80 cm s -1 respectively. Only very small fraction of these clumps survives the early stage of 
tidal destruction during the hierarchial clustering [3]- Nevertheless these survived clumps may 
provide the major contribution to the annihilation signal in the Galaxy [3 13 EH [231 123j- At a 
high redshift the SUSY neutralinos, considered as DM particles, may cause the efficient heating 
of the diffuse gas [25] due to annihilation in the dense clumps. 

One of the unresolved problem of DM clumps is a value of the central density or core radius. 
Numerical simulations give a nearly power density profile of DM clumps. Both the Navarro- 
Frenk-White (NFW) and Moore profiles give formally a divergent density in the clump center. 
A theoretical modeling of the clump formation ESEZlEij predicts a power-law profile of the 
internal density of clumps 

AntM = AAp()j) ", (1) 

where p and R are the mean internal density and a radius of clump, respectively, /3 — 1.8 — 2 
and = 0 at r > R. A near isothermal power-law profile (HD with f3 ~ 2 has been recently 

obtained in numerical simulations of small-scale clump formation {2Hj . In [3] the core radius 
x c — 0.3^ -2 has been obtained, where u is a relative height of the fluctuation density peak in 
units of dispersion at the time of energy-matter equality. This value is a result of the influence of 
tidal forces on the motion of DM particles in the clump at stage of formation. This estimate may 
be considered as an upper limit for the core radius or as the break-scale in the density profile, 
e. g., a characteristic scale in the Navarro-Frenk-White profile. It could be that a real core radius, 
where the density ceases to grow, is determined by the relaxation of small-scale perturbations 
inside the forming clump m • Another mechanism for core formation arises in the ‘meta-cold 
dark matter model’ due to late decay of cold thermal relics into lighter nonrelativistic particles 
with low phase-space density [Ml E2]. Here we consider the relative core radius x c = R c /R of 
DM clumps as a free parameter in the range 0.001 — 0.1. Correct approach includes a gradual 
mass loss of a systems [33l [Ml I36| . in particular, by small-scale DM clumps mm- 

2. Formation and Destruction of clumps in hierarchial clustering 

The process of hierarchical clustering and tidal destruction of DM clumps can be outlined in 
the following way. The DM clumps of minimal mass are formed first in the expanding Universe. 
The clumps of larger mass, which host the smaller ones are formed later, and so on. Some 
part of clumps are destroyed tidally in the gravitational field of their host clumps. In this 
Section we study the destruction of DM clumps in the process of hierarchical structuring long 
before the final galaxy formation. At small-mass scales the hierarchial clustering is a fast and 
rather complicated nonlinear process. A formation time of clump with an internal density p is 
t = (i^Peq/pY^teq, where k = 187r 2 and p eq = po(l + £eq) 3 is a cosmological density at the time 
of matter-radiation equality f e q; 1 + A?q = 2.35 x 10 4 Q m h 2 and pq = 1.9 x 10 _29 f i m h 2 g cm -3 . 
The index ‘eq’ here and throughout below refers to quantities at the time of matter-radiation 



equality t eq . The DM clumps of mass M can be formed from density fluctuations of different 
peak-height v = 5 eq /a eq (M), where a eq (M) is the fluctuation dispersion on a mass-scale M at 
the time t eq . A mean internal density of clump p is fixed at the time of clump formation and 
according to [38] equals p = np eq [ua eq (M )/ S c ] 3 , where 5 C = 3(127t) 2 / 3 /20 ~ 1.686. 

An internal energy of self-gravitating object increases in tidal interactions. This energy 
increase was calculated, e. g., in [3S] for the case of a star globular cluster in a spherical galaxy. 
By using the model of tidal heating from [39] . we determine now a survival time (or a time 
of tidal destruction) T of some chosen small-scale clump due to the tidal heating inside of a 
larger mass host clump. The motion of a clump would be rather complicated in the case of a 
fast hierarchical clustering of hosts. During a dynamical time in the host t^ yn ~ 0.5 (Gph) -1 ^ 2 , 
where ph is a mean internal density of the host, the chosen small-scale clump can belong to 
several successively destructed hosts. A clump trajectory in the host experiences successive 
turns accompanied by the “tidal shocks” [iCK 39] . For the considered small-scale clump with a 
mass M and radius R, the corresponding internal energy increase after a single tidal shock is 

Ajr 

A E~— 7l Gp h MR 2 , (2) 

where a numerical factor 71 ~ 1. Let us denote the number of tidal shocks per dynamical 
time fdyn by 72. A corresponding rate of clump internal energy growth is E = ^AE/td yn - A 
clump is destroyed in the host if its internal energy increase due to tidal shocks exceeds a total 
energy \E\ ~ GM 2 /2R. As a result, for a typical time T = T(p,ph ) of the tidal destruction of 
a small-scale clump with density p inside a more massive host with a density ph we obtain: 

T~\p,p h ) = E/\E\ ~ 4 7 l 72 G 1 / 2 / 0 fy i . (3) 

It turns out that a resulting mass function of small-scale clumps (see in this Section below) 
depends rather weakly on the value of 7172 . 

The probability of clump survival, determined as a fraction of the clumps with mass M 
surviving the tidal destruction in hierarchical clustering, is given by the exponential function 
e~ J with 

T \ ' Atfe 
~~ ^ 

Here At ^ is a difference of formation times t^ for two successive hosts, and summation goes over 
all chimps of intermediate mass-scales, which successively host the considered small-scale clump 
of a mass M. Changing the summation by integration in Q we obtain 


J(P, Pf) 


t f 


dt h Pi-Pf Pi t 2 

npiPh )- 1 p - 7 tr 


(5) 


where 

7 = 27i72K 1/2 G 1 ' /2 / 9g/ 2 t eq ^ 14(7 i7 2 /3), (6) 

and t, t±, tf, p, pi and pf are respectively the formation times and internal densities of the 
considered clump and of its first and final hosts. One may see from Eq. © that the first 
host provides a major contribution to the tidal destruction of the considered small-scale clump, 
especially if the first host density pi is close to p , and consequently e~ J <C 1 . 

Now we need to track the number of clumps M (originated from the density peak u) which 
enter some larger host during time intervals At± around each t± beginning from the time t of 





clump formation. A mass function of small-scale clumps (i. e., a differential mass fraction of DM 
in the form of clumps survived in hierarchical clustering) can be expressed as 
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(7) 


In this expression to is the Universe age and F(M , t) is a mass fraction of unconfined clumps (i. e., 
clumps, not belonging to more massive hosts) with a mass smaller than M at time t. According 
to [38], the mass fraction of unconfined clumps is F(M,t ) = erf (d c /[y/2a eq (M)D(t)\ j, where 
erf(x) is the error-function and D(t) is the growth factor normalized by D(t eq ) = 1. An upper 
limit of integration to i n Eh- m is not crucial and may be extrapolated to infinity because a 
main contribution to the tidal destruction of clumps is provided by the early formed hosts at 
the beginning of the hierarchical clustering. Two processes are responsible for time evolution of 
the fraction d 2 F/(dMdt ) for unconfined clumps in the mass interval dM: (i) the formation of 
new clumps and (ii) the capture of smaller clumps into the larger ones. Both these processes 
are equally efficient at the time when d 2 F/(dMdt ) = 0. Finally, we transform the distribution 
function (JTJ) to the following form: 


dM 


du 




v —v 2 /2 f / a ^ °eq (M) 

^ M7J dM 


dM, 


( 8 ) 


where 


2[r(l/3) - r(l/3,7)] 


(9) 


Here T(l/3) and r(l/3,7) are the Euler gamma-function and incomplete gamma-function, 
respectively. The function / 1 X 7 ) varies rather slowly in the interesting interval of 14 < 7 < 40, 
and one may use /i( 7 ) — 0.2 — 0.3. Physically the first factor v in (JHJ) corresponds to a more 
effective survival of high-density clumps (i.e., with large values of v) with respect to the low- 
density ones (with small values of v). Integrating Eq. (JHJ) over v, we obtain 


dM . .dM 

^ int —~ 0 . 02 (n + 3 )—. 


( 10 ) 


An effective power-law index n in Eq. m is determined as n = —3(1 + 2d log a eq (M)/d log M) 
and depends very weakly on M. 

In further calculations we use an interpolation fitting of the fluctuation dispersion cr eq (M) 
from m and PQ. The analysis of the WMAP data of the CMB anisotropy [42j reveals a 
power-law spectrum of initial perturbations with n p = 0.99 ± 0.04 in a good agreement with the 
canonical inflation value n p = 1.0. However, when the data from 2dF galaxy power-spectrum and 
Ly-a are included in the analysis, the best-fit favors in a softer spectrum with n p = 0.96 ± 0.02. 

The simple M -1 shape of the mass function GOD is in a very good agreement with the 
corresponding numerical simulations [29], but our normalization factor is a few times smaller. It 
must be stressed that a physical meaning of the survived clump distribution function £ dv dM/M 
is different from the similar one for the unconfined clumps, given by the Press-Schechter mass 
function dF/dM. 


3. Tidal destruction of clumps by the Galactic disk and stars 

Crossing the Galactic disc, a clump can be tidally destructed by the collective gravitational field 
of stars in the disc. This phenomenon is similar to the destruction of a globular cluster by the 









“tidal shocking” in the Galactic disc [23j. The kinetic energy gain of a DM particle with respect 
to the center of clump after one crossing of the Galactic disk is @3] 
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( 11 ) 


where m is a constituent DM particle mass, Az is a vertical distance (orthogonal to the disk 
plane) of a DM particle with respect to the center of clump, v Z}C is a vertical velocity of clump 
with respect to the disk plane at the moment of disk crossing and A{a) is the adiabatic correction 
factor. A gravitational acceleration near the disk plane is g m {r ) = 2irGcr s {r), where we use an 
exponential model for a surface density of disk. The factor A{a) in (HID describes the adiabatic 
protection from slow tidal effects [43]. In [33] the following fitting formula was proposed: 
A{a) = (1 + a 2 ) _3//2 . Here the adiabatic parameter a = uiTd, where u> is an orbital frequency 
of DM particle in the clump, ~ Hd/v ZyC is an effective duration of gravitational tidal shock 
produced by the disk with a half-thickness Hd- 

As a representative example we consider the isothermal internal density profile of DM clump 


Pint (0 


1 


v 2 

Got 


4-7T Gr 2 


( 12 ) 


with a cutoff at the virial radius R: p(r ) = 0 at r > R. A corresponding mass profile of 
clump is M(r ) = AR{r/R), where Mi is an initial mass of clump at the epoch of Galaxy 
formation. With this mass distribution a circular velocity inside a clump is independent of 
radius, n rot = (GM(r)/r)^ 2 = {GAR/R ) 1 / 2 . A gravitational potential corresponding to the 
density profile CED is <f>(r) = v 2 ot [log(r/i?) — 1]. Let us define a dimensionless energy of the 
DM particle e = E/{mv 2 ot ) and gravitational potential i/{r) = 4>{r)/v/ ot = In {r/R) — 1. An 
internal density profile Pint( r ) and the distribution function of DM particles in the clump f c i(e) 
are related by the integral relation [45]] 


o _ 

Pint(r) = 2 5/2 vr J \J£ — V’( r ) /ci(^) de, (13) 

4>(r) 

The corresponding isothermal distribution function is f c i(e) oc exp(—2e). 

By using the hypothesis of a tidal stripping of outer layers of a DM clump, we see that a 
tidal energy gain de causes the stripping of particles with energies in the range — Se < e < 0. A 
corresponding variation of density at radius r is 

o _ 

Sp(r) = 2 5/2 7T J yje - ip{r) / c i(e) de. (14) 

—Se 


In this equation the tidal energy gain CD by different DM particles is averaged over angles, so 
as ((A^) 2 ) = r 2 /3. A resulting total mass loss by DM clump during one crossing of the Galactic 
disk is 


8M = —47r / r 2 5p{r ) dr. 

Jo 


(15) 


Let us specify the dimensionless quantities 
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where pi = 3 Mj/( 47 ri? 3 ) is a initial mean density of clump. For the most parts of clumps Qd <C 1 
with a typical value Qd ~ 0.03. The fitting formula for the mass loss of clump during one passage 
through the Galactic disk is 


(jr) d ~ -°- 13Qdexp (~ L58 ^ /2 ) • (i7) 

Let us choose some particular clump moving in the spherical halo with an orbital “inclination” 
angle 7 between the normal vectors of the disk plane and orbit plane. The orbit angular velocity 
at a distance r from the Galactic center is dcp/dt = J/(mr 2 ), where J is an orbital angular 
momentum of a clump. A vertical velocity of a clump crossing the disk is 


J . 

v z ,c = -sin 7 , 

mr s 


(18) 


where r s is a radial distance of crossing point from the Galaxy center. There are two crossing 
points (with different values of r s ) during an orbital period. 

The standard Navarro-Frenk-White profile of the DM Galactic halo is 


Pair) 


Po 
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(19) 


where L 


45 kpc, po = 5 x 10 6 Mq kpc 3 . 


It useful to introduce the dimensionless variables: 
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where 4>o = ~^GpoL 2 , E OT b is a total orbital energy of a clump. With these variables the 
density profile of the halo (fT9l) is written as pu{x) = [x(l + x) 2 ] -1 . A gravitational potential 
tp{x), corresponding to density profile (fl9ll is ip{x) = 1 — log(l + x)/x. The relation between 
density profile Ph{x) and the distribution function is given by the same equation (1131) with an 
obvious substitution f c \ => F(e), where the distribution function F{e) for a halo profile (fT9l) can 
be fitted as in [45]. 

Let us denote the interval of time for motion from x m i n to x max and back as T c (x,e,p), and 
the angle of orbital precession during the time TJ 2 as cp. The orbital period is longer than T c 

and is given by T< = T c ^1 + (p/irj . Choosing a time interval AT much longer than a clump 
orbital period T t . but much shorter than the age of the Galaxy to, i. e., T t <C AT <C to, we dehne 
an averaged rate of mass loss by a selected clump under influence of tidal shocks in successive 
disk crossings: 

M \ dt ) d ~ AT ) d ( j 

where ( 5M/M)d is given by ()17h and summation goes over all successive crossing points (odd 
and even) of a clump orbit with the Galactic disk during the time interval AT. The values of 
the g m and v Zfi both depend on the radius x = r/L. One simplification in calculation of (12TT) 
follows from the fact that a velocity of orbit precession is constant. For this reason the points of 
successive odd crossings are separated by the same angles 4>. The same is also true for successive 
even crossings. Using this simplification we transform the summation in (1211) to integration: 
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where d(j>/dx is an equation for clump orbit in the halo. 

During a single close encounter of a DM clump with a star, the energy gain of a constituent 
DM particle in the clump with respect to clump center is [6]: 
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( 22 ) 


where m* is a star mass, l is an impact parameter, i; re i is a relative star velocity with respect 
to a clump, A z = r cos ip, r is a radial distance of a DM particle from the clump center and -0 
is an angle between the directions from the clump center to the DM particle and to the point 
of closest approach of a star. Using the same method as in the Section [3] we calculate a relative 
mass loss by clump ( SM/M) s during a single encounter with a star and obtain the same fitting 
formula as (ED but with a substituting the dimensionless parameters, Qd => Qs and Sd =>- S s , 
where 
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(23) 


where t s — l/v re i. 

A DM clump acquires the maximum energy gain during a single encounter with a star when 
impact parameter l ~ R. Integrating over all impact parameters l > R, we calculate an averaged 
rate of mass loss by a clump during successive encounters with stars: 


OO Xmax 

J_ (M£\ ~ 1 [ 2ldd i f dsn*{s)v rel / 6M_\ 
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(24) 


where n*(r) is a radial number density distribution of stars in the bulge and halo. 


4. Surviving fraction of clumps 

From Eq. ED it is seen that the tidal forces influence mainly the outer part of clump (where A z 
is rather large). Further we use our basic assumption that only outer layers of a clump undergo 
the tidal stripping, while the inner parts of a clump are unaffected by tidal forces. Thus we 
assume that a clump mass M = M(t) and radius R = R(t) are both diminishing in time due to 
the tidal stripping of outer layers but its internal density profile remains the same as given by 
Eq. (1121b e. g., for the isothermal density profile M(t) oc R(t) and p(t) oc M(t)~ 2 . Combining 
together the rates of mass loss (12TT) and (|24l) due to the tidal stripping of a clump by the disk 
and stars respectively, we obtain the evolution equation for a clump mass: 


dM _ fdM\ / dM\ 
dt V dt ) d + V dt ) s 


(25) 


In the following we solve this equation numerically starting from the time of Galaxy formation 
at to ~ to up to the present moment to- 

The most important astrophysical manifestation of DM clumps is a possible annihilation 
of the constituent SUSY DM particles. The crucial point is a dominance of the central core 
of a clump in annihilation signal if chimps have a steep enough density profile. Namely, 
annihilation of DM particles in a clump core prevail in a total annihilation rate in a single 
clump with a power-law density profile CD if /3 > 3/2 and x c = R c /R 1. More specifically, 
the quantity N oc J)' o iirr 12 dr'pf nt (r') practically does not depend on r, if r 3> tq. As a result the 
annihilation luminosity of DM clump with approximately isothermal density profile (/? ~ 2) will 
be nearly constant under influence of tidal stripping until a clump radius diminishes to its core 
radius. In other words, in the nowadays Galaxy the remnants of tidally stripped chimps with 








x c < //(to) -C 1, where p(t) = M(f)/Mj and to — 10 10 yrs is the Galaxy age, obeys the evolution 
equation (1251) and have the same annihilation luminosity as their progenitors with /j = 1. 

By using evolution equation (1251) we now calculate the probability P of the survival of clump 
remnant during the lifetime of the Galaxy. Let us choose some arbitrary point in the halo with 
a radius-vector r and an angle a with a polar axis of the Galactic disk. Only the clump orbits 
with inclination angle 7 t/2 — a < 7 < it/2 pass through this point. A survival probability for 
clumps can be written now in the following form 

1 sin a 1 

P(x, a) = - f dp f dcos^ f de [e — ^(x)] 1 ^ 2 F(e)0[n(to) — x c ], (26) 

c)\ocj sin O! 1 / ,/ j 

0 0 ^;(x) 

In this equation p(x) is a density profile of the halo from (1 19 1) . p = cos 0, 6 is an angle between 
the radius-vector r and the orbital velocity of clump, 0 is the Heaviside function, i/>(x) is the 
halo gravitational potential from (1191) . F(e) is a distribution function of clumps in the halo, //(to) 
depends on all variables of the integration and x c = R c /R is an initial value of clump core. The 
function //(to) is calculated from numerical solution of evolution equation ([25]) . If //(to) > x c , 
the clump remnant is survived through the tidal destruction by both the disk and stars. The 
annihilation rate in this remnant would be the same as in the initial clump. On the contrary in 
the opposite case, when //(t^) < x c , the clump is totally destructed because (i) the core is not a 
dynamically separated system and composed of particles with extended orbits, and because (ii) 
a nearly homogeneous core is destructed easier than a similar object with the same mass but 
with a near isothermal density profile. 

According to theoretical model f[3) and numerical simulations [29], a differential number 
density of small-scale clumps in the co-moving frame in the Universe is n(M)dM oc dM/M 2 . 
This distribution is calculated in [ 6 ]. The damping of small-scale perturbations with M < M m ; n 
provides an additional factor exp[— (M/M m i n ) 2 / 3 ] responsible for the fading of distribution 
at small M. The result of numerical simulations |29| can be expressed in the form of a 
differential mass fraction of the DM clumps in the Galactic halo f(M)dM ~ n(dM/M), where 
k — 8.3 x 10~ 3 . The analytical estimation ( 1101 ) gives approximately s~4x 1CD 3 for the mass 
interval 10~ 6 Mq < M < IMq. The discrepancy by the factor ~ 2 may be attributed to the 
approximate nature of our approach as well as to the well known additional factor 2 in the 
original Press-Schechter derivation of the mass function. In the latter case one must simply 
multiply equation (0 by factor 2 . 

5. Amplification of annihilation signal 

A local annihilation rate is proportional to the square of DM particle number density. A number 
density of DM particles in clump is much large than a corresponding number density of the diffuse 
(not clumped) component of DM. For this reason an annihilation signal from even a small fraction 
of DM clumps can dominate over an annihilation signal from the diffuse component of DM in 
the halo. We consider here the Harrison-Zeldovich initial perturbation spectrum with power 
index n p = 1 as a representative example. The value of n p is not exactly fixed by the current 
observations of CMB anisotropy. In the case of n p < 1, the DM clumps are less dense, and a 
corresponding amplification of annihilation signal would be rather small [3]. 

The gamma-ray flux from annihilation of diffuse distribution m of DM in the halo is 
proportional to 

^max(C) 

I U = I Ph(£) dx, (27) 

0 

where the integration is over r goes along the line of sight, £((j r ) = (r 2 + Tq — 2 rr® cos C ) 1//2 
is the distance to the Galactic center, r max (£) = ( R ^ — 7 q sin 2 (j) 1 / 2 + r Q cos() is a distance to 



the external halo border, £ is an angle between the line of observation and the direction to the 
Galactic center, Rn is a virial radius of the Galactic halo, r© = 8.5 kpc is the distance between 
the Sun and Galactic center. The corresponding signal from annihilations of DM in clumps is 
proportional to the quantity [3] 

^max(C) 

Ici = S J dx j f(M)dMpp H (Z)P(S,p), (28) 

0 Mmin 

where p(M ) is the mean density of clump. The function S depends on the clump density profile 
and core radius of clump |3] and we use S — 14.5 as a representative example. The observed 
amplification of the annihilation signal is defined as 77 (C) = (4 c i+/h)/^h is calculated numerically 
in Ei for the case x c = 0.1. It tends to unity at £ —» 0 because of the divergent form of the 
halo profile () 19 1) . The annihilation of diffuse DM prevails over signal from chimps at the the 
Galactic center. The observed signal is obtained by integration along the line of sight and the 
effect of clumps destruction at the Galactic center is masked by the signal from another regions 
of the halo. This amplification of an annihilation signal is called a “boost-factor” m- 

6. Anisotropy of the observed signal 

The usual assumption in calculations of DM annihilation is a spherical symmetry of the Galactic 
halo. In this case an anisotropy of annihilation gamma-radiation is only due to off-center position 
of the Sun in the Galaxy. Nevertheless, a principal significance of the halo nonsphericity for the 
observed annihilation signal was demonstrated in |48j . According to observations, the axes of the 
Galactic halo ellipsoid differ most probably no more than 10 — 20%, but even a much more larger 
difference of axes, up to a factor 2, can not be excluded @9, 50]. This leads to more than an 
order of magnitude uncertainty in the predicted annihilation flux from the Galactic anti-center 
direction [T8J. It was shown ee that (i) small-scale DM clumps dominate in the generation 
of annihilation signal and (ii) the Galactic stellar disk provides the main contribution to the 
tidal destruction of chimps at r > 3 kpc, i. e., outside the central bulge region. The detectors 
at the Fermi satellite are sensitive to anisotropy up to 0.1% level m This provide a hope to 
discriminate the anisotropic DM annihilation signal from the diffuse gamma-ray backgrounds. 
To calculate the anisotropy we used a simplified criterium for a tidal destruction of clump: clump 
is destructed if a total tidal energy gain XX A E)j after several tidal shocking becomes of order 
of initial binding energy of a clump \E\, i. e., A E = XX A E)j ~ \E\, where summation goes over 

j 

the successive disk crossings (or encounters with stars). In this model the survival probability 
for clumps can be written now in the following form 

1 sin a 1 

P(r, a) = ——-— / dp [ dcosy f de [e — iplx)] 1 ^ 2 F(e) ( 29 ) 

p(x) sm a J J J 

0 0 V'M 

where the definitions are the same as in (1261) except for exponential factor for clumps destruction 
e —AE/\E\' j n gjj thg triple integral from (|29l) for survival probability P(r,a ) was numerically 
calculated. As a representative example we consider the Earth-mass clumps M = 10~ 6 Mq 
originated from 2cr density peaks in the case of power-law index of primordial spectrum of 
perturbations n p = 1. The mean internal density of these clumps is p ~ 7 x 10 -23 g cm -3 . The 
annihilation signal for the Galactic disk plane and for the orthogonal vertical plane (passing 
through the Galactic center) as a function of angle £ between the observation direction and the 
direction to the Galactic center was calculated numerically in mi- The difference of the signals 
in two orthogonal planes at the same ( can be considered as an anisotropy measure. Defined as 
5 = (I 2 — Ii)/Ii, it has a maximum value 6 — 0.09 at £ ~ 39°. 



7. Conclusion 

We demonstrate that the cores of the small-scale DM clumps (or clump remnants) survive in 
general during clump formation and their tidal destruction by stars in the Galaxy. These small- 
scale DM clumps may provide the major contribution to the annihilation signal (in comparison 
with the diffuse DM) in the Galactic halo. The amplification (boost-factor) can reach 10 2 or 
even 10 3 depending on the initial perturbation spectrum and minimum mass of clumps. 
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